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Superfluidity of A and I]~ admixed in neutron star (NS) cores is investigated realistically 
for hyperon (Y')-mixed NS models obtained using a G-matrix-based effective interaction ap- 
I proach. Numerical results for the equation of state (EOS) with the mixing ratios of the 

04 , respective components and the hyperon energy gaps including the temperature dependence 

are presented. These are meant to serve as physical inputs for y-cooling calculations of NSs. 
CN ■ By paying attention to the uncertainties of the EOS and the YY interactions, it is shown 

^ I that both A and S~ are superfluid as soon as they appear although the magnitude of the 

CO , critical temperature and the density region where superfluidity exists depend considerably 

\l ' on the YY pairing potential. Considering momentum triangle condition and the occurrence 

' of superfluidity, it is found that a so-called "hyperon cooling" (neutrino-emission from direct 

, Urea process including Y) combined with F-superfluidity may be able to account for obser- 

^— »' . vations of the colder class of NSs. It is remarked that /1-hyperons play a decisive role in the 

hyperon cooling scenario. Some comments are given regarding the consequences of the less 
>— ^ ' attractive AA interaction recently suggested by the "NAGARA event" ^yiHe. 

§1. Introduction 

O 

pH I Neutron stars (NSs) are primarily composed of neutrons, with a small amount 

of protons, together with electrons and muons to assure the charge neutrality of the 
system. With the increase of the total baryon density p toward the central region, 
r> I however, hyperons (Y), such as A, E~ and S~ , begin to appear as new components. 

■ This is because the chemical potential of the predominant neutron component gets 

higher with p, and eventually it becomes energetically economical for the system to 
replace neutrons at the Fermi surface by Y through a strangeness non-conserving 
weak interaction, in spite of the cost of the higher rest-mass energy. The population 
of Y increases with increasing p and hyperons, in addition to nucleons, become 
important constituents of NS cores. 

Internal composition and structure strongly affect the thermal evolution of NSs, 
and this problem has been extensively studied from various points of view.^^^^*^^ 
The participation of Y provides a much faster cooling mechanism due to the very 
efficient neutrino-emission processes, ^^-^ the so-called "hyperon direct Urea" (abbre- 
viated as y-DUrca; e.g., A ^ p-\-l-\-i>i^ E~ A+l + i>i and their inverse processes; 
Z = or /x^), compared to the standard cooling mechanism, called the "modified 
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Urea" (MUrca; e.g., n + n n + p + I + ui, n + p -^p + p+ l + ui and their 
inverse processes). The z^-emission rate for hyperon coohng (y-DUrca) is larger by 
5-6 orders of magnitude than that for standard coohng (MUrca), because in the 
latter an additional nucleon has to join in the /3-decay process in order to satisfy 
momentum conservation for the reactions. This rapid hyperon cooling is of partic- 
ular interest in relation to the surface temperature observations of NSs. That is, 
recent observations suggest that there are at least two classes of NSs, hotter ones 
and colder ones, and these findings have motivated studies on cooling scenarios to 
explain the two classes consistently. Standard cooling, together with a frictional 
heating mechanism, can explain the hotter ones, but not the colder ones. Then 
the fast non-standard hyperon cooling shows up as a candidate to explain the colder 
ones. The direct action of hyperon cooling, however, leads to a serious problem of 
"too rapid cooling" incompatible with observations and therefore some suppression 
mechanism for an enhanced z^-emission rate is necessary. The most natural choice 
as a possible mechanism to play this suppressive role is hyperon superfluidity. 

In this paper, we study whether hyperons admixed can be superfluids and show 
that indeed they can.^)"^) This implies that there is strong possibility that the hy- 
peron cooling scenario, combined with the y-superfiuidity, is the mechanism re- 
sponsible for fast non-standard cooling compatible with the data of the colder class 
of NSs. In our preceding works, ^'^^^^^^ we studied the hyperon energy gap Ay re- 
sponsible for y-superfluids using a simplified treatment in which the p-dependent Y- 
mixing ratios yy(= Py / p) relevant to the Fermi momentum qpy and the /9-dependent 
effective-mass parameter my are taken as parameters. In the present work, we use 
/9-dependent yvip) and my(/3) based on the Y -mixed NS models in order to derive 
a realistic Ay in a p-dependent way. Also, we present profile functions represent- 
ing the temperature-dependence of Ay, which is necessary for treating the thermal 
evolution of NSs, and attempt to provide the basic physical quantities for cooling 
calculations of a realistic level. 

This paper is organized as follows. In §2, our NS models with F-mixing, calcu- 
lated using a G-matrix effective interaction approach, are presented, together with 
the composition, the equation of state (EOS) and the physical quantities neces- 
sary for the energy gap calculations. The energy gap equation including a finite- 
temperature (T > 0) case is treated in §3, where some comments are given on the 
impotant ingredients in the gap equation, i.e., my(p) and the YY pairing inter- 
actions in comparison with the case of nucleons. Results for Ay and the profile 
functions are given in §4, together with discussion regarding the momentum triangle 
condition for F-DUrca processes and the superfluid suppression of z/-emissivities. In 
addition, in §5, we discuss how a less attractive AA interaction inferred from ^yjHe 
observed recently affects the realization of F-superfluidities. The last section, §6, 
contains a summary and remarks. 

§2. Neutron Star Models with Hyperons 

In this section, we calculate the EOS of Y-mixed NS matter responsible for Y- 
mixed NS models and also obtain the p-dependent yy and my of hyperons necessary 
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for the energy gap calculations. For simplicity, we ignore the mixing of the S~ com- 
ponent and restrict our investigation to Y = {A, U~} as the hyperon components, 
since it would be less likely for S" to appear due to the larger mass {m^- = 1321 
MeV compared to my[ = 1116 MeV and m^- = 1192 McV). Wc treat Y-mixed NS 
matter composed of n, p, A, U~ , e~ and using a G-matrix based effective interac- 
tion approach, with attention to the choice of YN and YY interactions compatible 
with hypernuclear data:^^)~^^^ 

1. Wc construct the effective YN and YY local potentials, Vy^ and Vyy, in a p- 
and ?/y-dependent way, from the G- matrix calculations performed for {n + Y} 
matter with a given p and yy- In these calculations, we use the YN and YY 
interactions from the Nijmegen model D hard core potential^e) (NHC-D) with 
a slight modification in the S'-statc AN part (NHC-Dm) as a reasonable choice 
(details are referred to Ref. 23)). From a view of A effective-mass parameter 
{m\) controlling the realization of yl-superfiuidity, which is discussed in the next 
section, the choice of this NHC-Dm is appropriate since it reproduces ~ 0.8 
infered from hypernuclear data.^^^ 

2. For the effective NN local potential Vnn, we use Vrsc supplemented by VtnIj 
i.e., Vnn = Vrsc + Vtni where Vrsc^^^ is the effective two-nucleon poten- 
tial constructed from the G-matrix calculations in asymmetric nuclear mat- 
ter with the Reid soft-core potential and Vtni^^^ is a phenomenological three- 
nucleon interaction (TNI) of the Lagaris-Pandharipande type.^''^'^^^ The po- 
tential Vtni(= ^tna + Vtnr) consists of two parts, the attractive Vtna and 
the repulsive VtnRj and is expressed in the form of a two-body potential with 
a p-dependence. At high densities the contribution from Vtna is minor, while 
that from Vtnr dominates, and TNI brings a strong repulsion with increasing 
p. The parameters inherent in Vtni are determined so that Vtni can reproduce 
the empirical saturation properties of symmetric nuclear matter (binding en- 
ergy Eb = —16 MeV and saturation density po = 0.17 nucleons/fm^) and the 
nuclear incompressibility k. 

3. On the basis of Vnn, Vyn and Vyy, we calculate the fractions yi for the re- 
spective components (i = n, p, A, E~ , e~ and fi^) in /3-equilibrium, under 
the conditions of charge neutrality, chemical equilibrium and baryon number 
conservation, and derive the EOS and my(p) of F-mixed NS matter. 

Using the EOS thus obtained, we solve the so-called TOV equation and obtain 
NS models with a y-mixed core. In the calculations, extended use of our G-matrix- 
based effective interactions to high-density regime is inevitable, and this introduces 
uncertainties with increasing p. These uncertainties, however, would not be serious 
because the p-dependence of E (the energy per particle), namely the EOS, is similar 
to that obtained with more realistic treatments, as shown in Fig.l for the TNI6 case, 
for example. Our NS models are specified by the parameter k relevant to Vnn, which 
is a measure of the stiffness of the nuclear-part EOS (A'^-part EOS). We consider three 
cases; k = 250 MeV (hereafter, referred to as TNI2), 300 MeV (TNI3) and 280 MeV 
(TNI6). TNI2 is a soft EOS case, TNI3 is a stiffer EOS case and TN16 is between the 
two, giving respectively the maximum mass M^aax of NSs sustained by the EOS as 
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Afmax ^ 1.62Mq, I.88M0 and l.TSM© for normal NSs without F-mixing. In Fig.l, 
we compare our EOSs for neutron matter with those by the Illinois group currently 
regarded as most realistic. We stress that our EOSs cover almost the entire range 
from soft to stiff cases allowed by realistic EOSs derived using a more sophisticated 
many-body approach. ^^^'^^^ 



TNI3 UV14+TNI 




Fig. 1. Energy per particle E versus density 

p, corresponding to the EOS of neutron 
matter (po = 0.17 nucleons/fm^ being 
the nuclear density). The solid curves 
labeled TNI2, TNI6 and TNI3 are re- 
spectively soft, intermediate and the stiff 
EOSs considered here. The EOSs from 
the Illinois group, i.e., AV18^^' (soft- 
est), UV14+TNI^^^ (intermediate) and 
AV18+(5v+UIX-''-''^ (stiffest), are also dis- 
played by the dotted lines. It is seen that 
our EOSs are well within the range from 
the softest to stiffest EOSs currently taken 
as realistic. 
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Fig. 2. Mass M versus central density pc for 
neutron stars (NSs) corresponding to the 
EOS of hyperon (y)-mixed neutron star 
matter (TNI2u, TNI6u and TNI3u), where 
repulsion from the three-nucleon interac- 
tion (TNI) is introduced universaly (see the 
main text for details). The crosses denote 
the maximum mass points. The threshold 
density pt{Y) of V-mixing and the density 
region, in which there exists yl-superfluid 
for ND-Soft pairing interaction, are also 
shown in order to indicate the M and EOS 
dependences of the yl-superfluid phase. 



Table I. Profile (mass Mmax, radius R and central density pc) of the maximum-mass NSs depicted 
in Fig.2 for EOSs from soft (TNI2u) to stiff (TNI3u). The thresold density pt{Y) indicating the 
portion of F-mixed core is also listed for the A and E~ cases. 



EOS 


Pt{A)/po 


PtiS-)/po 


Mmax /Mo 


R/km 


Pc/po 


TNI2u 


4.01 


4.06 


1.52 


8.43 


11.08 


TNI6u 


4.02 


4.06 


1.71 


9.16 


9.07 


TNISu 


4.01 


4.01 


1.83 


9.55 


8.26 
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Table II. (a) EOS of Y- mixed NS matter with mixing ratio yi of the constituents (i = n, p, A, , 
e~ , iJb~) for the case of TNI2u. The energy density e is in units of MeV/fm^ and the pressure 
P is in units of dyn/cm^. 



p/po 


Vn 


Vp 


VA 


Ve- 




y^- 


e 


P 


0.5 


0.9684 


0.0316 


0. 


0. 


0.0316 


0. 


80.5 


5.51E+32 


1.0 


0.9499 


0.0501 


0. 


0. 


0.0451 


0.0050 


161.9 


3.67E+33 


1.5 


0.9409 


0.0591 


0. 


0. 


0.0464 


0.0128 


244.9 


1.13E+34 


2.0 


0.9380 


0.0620 


0. 


0. 


0.0453 


0.0168 


330.1 


2.40E+34 


2.5 


0.9371 


0.0629 


0. 


0. 


0.0439 


0.0191 


417.6 


4.23E+34 


3.0 


0.9368 


0.0632 


0. 


0. 


0.0426 


0.0206 


507.7 


6.70E+34 


3.5 


0.9368 


0.0632 


0. 


0. 


0.0415 


0.0216 


600.9 


9.92E+34 


4.0 


0.9371 


0.0629 


0. 


0. 


0.0406 


0.0224 


697.2 


1.40E+35 


4.5 


0.8561 


0.0862 


0.0262 


0.0315 


0.0354 


0.0193 


796.8 


1.80E+35 


5.0 


0.7633 


0.1127 


0.0574 


0.0667 


0.0301 


0.0159 


899.2 


2.25E+35 


5.5 


0.6824 


0.1354 


0.0859 


0.0963 


0.0258 


0.0133 


1004.8 


2.78E+35 


6.0 


0.6135 


0.1545 


0.1109 


0.1211 


0.0223 


0.0111 


1113.5 


3.38E+35 


6.5 


0.5556 


0.1701 


0.1327 


0.1416 


0.0193 


0.0092 


1225.6 


4.08E+35 


7.0 


0.5071 


0.1827 


0.1518 


0.1584 


0.0167 


0.0076 


1341.3 


4.87E+35 


7.5 


0.4660 


0.1929 


0.1687 


0.1723 


0.0144 


0.0062 


1460.7 


5.77E+35 


8.0 


0.4308 


0.2013 


0.1838 


0.1840 


0.0123 


0.0049 


1584.1 


6.77E+35 


9.0 


0.3734 


0.2143 


0.2096 


0.2027 


0.0088 


0.0028 


1843.3 


9.10E+35 


10.0 


0.3277 


0.2243 


0.2309 


0.2170 


0.0060 


0.0013 


2120.3 


1.19E+36 


11.0 


0.2898 


0.2329 


0.2485 


0.2288 


0.0038 


0.0003 


2416.3 


1.52E+36 


12.0 


0.2570 


0.2410 


0.2631 


0.2388 


0.0022 


0.0000 


2732.5 


1.90E+36 


13.0 


0.2272 


0.2492 


0.2754 


0.2481 


0.0011 


0.0000 


3069.9 


2.34E+36 


14.0 


0.1983 


0.2581 


0.2859 


0.2577 


0.0004 


0.0000 


3429.3 


2.82E+36 


15.0 


0.1675 


0.2689 


0.2948 


0.2688 


0.0001 


0.0000 


3811.6 


3.36E+36 


(b) Same as (a), but for the TNI6u case. 


p/po 


Vn 


Vp 




Vs- 


y^- 




e 


P 


0.5 


0.9682 


0.0318 


0. 


0. 


0.0318 


0. 


80.5 


5.87E+32 


1.0 


0.9473 


0.0527 


0. 


0. 


0.0468 


0.0058 


162.0 


4.00E+33 


1.5 


0.9369 


0.0631 


0. 


0. 


0.0488 


0.0143 


245.4 


1.31E+34 


2.0 


0.9341 


0.0659 


0. 


0. 


0.0475 


0.0184 


331.2 


2.91E+34 


2.5 


0.9341 


0.0659 


0. 


0. 


0.0455 


0.0203 


420.2 


5.26E+34 


3.0 


0.9349 


0.0651 


0. 


0. 


0.0437 


0.0214 


512.6 


8.47E+34 


3.5 


0.9358 


0.0643 


0. 


0. 


0.0421 


0.0221 


608.8 


1.27E+35 


4.0 


0.9366 


0.0634 


0. 


0. 


0.0408 


0.0226 


709.2 


1.80E+35 


4.5 


0.8550 


0.0869 


0.0261 


0.0321 


0.0354 


0.0193 


814.0 


2.37E+35 


5.0 


0.7609 


0.1136 


0.0578 


0.0677 


0.0300 


0.0159 


923.0 


3.02E+35 


5.5 


0.6797 


0.1363 


0.0865 


0.0975 


0.0257 


0.0132 


1036.4 


3.79E+35 


6.0 


0.6110 


0.1553 


0.1116 


0.1222 


0.0221 


0.0109 


1154.5 


4.68E+35 


6.5 


0.5534 


0.1707 


0.1334 


0.1425 


0.0191 


0.0091 


1277.6 


5.71E+35 


7.0 


0.5052 


0.1831 


0.1525 


0.1591 


0.0165 


0.0075 


1406.1 


6.89E+35 


7.5 


0.4645 


0.1932 


0.1694 


0.1729 


0.0142 


0.0061 


1540.1 


8.23E+35 


8.0 


0.4296 


0.2015 


0.1845 


0.1845 


0.0122 


0.0048 


1680.0 


9.73E+35 


9.0 


0.3725 


0.2144 


0.2102 


0.2029 


0.0087 


0.0028 


1978.6 


1.32E+36 


10.0 


0.3271 


0.2244 


0.2313 


0.2172 


0.0059 


0.0013 


2304.1 


1.75E+36 


11.0 


0.2893 


0.2330 


0.2489 


0.2289 


0.0038 


0.0003 


2658.3 


2.25E+36 


12.0 


0.2566 


0.2411 


0.2635 


0.2389 


0.0022 


0.0000 


3043.3 


2.83E+36 
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(c) Same as (a), but for the TNI3u case. 



p/Po 


Vn 


Vp 




Vs- 


Vc- 




e 


P 


0.5 


0.9681 


0.0319 


0. 


0. 


0.0319 


0. 


80.5 


6.17E+32 


1.0 


0.9454 


0.0546 


0. 


0. 


0.0481 


0.0065 


162.1 


4.21E+33 


1.5 


0.9334 


0.0666 


0. 


0. 


0.0509 


0.0158 


245.6 


1.42E+34 


2.0 


0.9302 


0.0698 


0. 


0. 


0.0498 


0.0201 


332.0 


3.26E+34 


2.5 


0.9306 


0.0694 


0. 


0. 


0.0475 


0.0219 


422.0 


6.01E+34 


3.0 


0.9321 


0.0679 


0. 


0. 


0.0452 


0.0226 


516.0 


9.80E+34 


3.5 


0.9339 


0.0661 


0. 


0. 


0.0432 


0.0230 


614.5 


1.48E+35 


4.0 


0.9354 


0.0646 


0. 


0. 


0.0415 


0.0231 


718.1 


2.11E+35 


4.5 


0.8457 


0.0907 


0.0274 


0.0362 


0.0353 


0.0192 


826.8 


2.79E+35 


5.0 


0.7516 


0.1172 


0.0594 


0.0718 


0.0297 


0.0157 


940.6 


3.59E+35 


5.5 


0.6712 


0.1395 


0.0881 


0.1012 


0.0254 


0.0129 


1059.8 


4.54E+35 


6.0 


0.6038 


0.1578 


0.1132 


0.1253 


0.0218 


0.0107 


1184.9 


5.65E+35 


6.5 


0.5475 


0.1726 


0.1349 


0.1449 


0.0189 


0.0089 


1316.2 


6.94E+35 


7.0 


0.5004 


0.1846 


0.1539 


0.1611 


0.0163 


0.0073 


1454.2 


8.41E+35 


7.5 


0.4605 


0.1943 


0.1707 


0.1744 


0.0140 


0.0059 


1599.1 


l.OlE+36 


8.0 


0.4263 


0.2023 


0.1857 


0.1857 


0.0120 


0.0047 


1751.3 


1.19E+36 


9.0 


0.3702 


0.2148 


0.2113 


0.2037 


0.0085 


0.0027 


2079.3 


1.63E+36 


10.0 


0.3255 


0.2246 


0.2322 


0.2177 


0.0058 


0.0012 


2440.9 


2.17E+36 



As a peculiar property of the y-mixed NS models, we remark that the y-mixing 
causes a dramatic softening effect on the EOS and consequently Mmax for y-mixed 
NSs is greatly reduced. For example, Mmax — 1.62M0 for the case without 
y is reduced to M ~ l.OSM© for the case with Y , contradicting the condition that 
Mmax should be larger than Mobs (PSR1913+16)=1.44Mq, the observed NS mass for 
PSR1913+16. This inconsistency between theory and observation cannot be resolved 
by using a stronger A^A^ repulsion, as seen from Afmax — l.OSM© — ^ l.O9Af0 — ^ 
I.IOMq for TNI2^TNI6 ^TNI3 with increasing k. This is because as the A-part 
EOS becomes stiffer, the y-mixed phase develops from lower densities, and as a 
result, the softening effect becomes stronger, making the enhanced A^A^ repulsion 
ineffective. The problem that the relation Mmax < Mobs is predicted for y-mixed 
NSs is also seen in other approaches, ^^^'^^^ and it is encountered in an almost model- 
independent way, which interestingly suggests that some "extra repulsion" must exist 
in hypernuclear systems, namely, yA" and YY interaction parts. As one possible such 
repulsion, we try to introduce a repulsion from the three-body force Vtnr in TNI also 
into the yA and YY parts, as well as the AA^ part (hereafter called the "universal 
inclusion of TNI" and denoted simply by TNIu),^^^"^^^ because the importance of 
the three-body interaction is well recognized for nuclear systems and should not be 
restricted to nuclear systems (A"A^ part). When this is done, we have a moderate 
softening, not a dramatic one, as before, and we obtain a larger Mmax that satisfies 
the condition M^ax > 1.44M0; M^ax ^ 1.52M0, I.TIM© and 1.83Mq for TNI2u, 
TNI6u and TNI3u, respectively. We also note that the realization of a y-mixed 
phase is pushed to the higher density side; the threshold density ptiY) ~ (2.5 — 3)po 
for no extra repulsion whereas ptiY) ~ 4po for the universal inclusion of the three- 
body repulsion. These results suggest that ptiX) is not so low as currently believed 
{pt{Y) ~ 2pq) but rather high (~ 4/9o), as far as consistency between the y-mixed 
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EOS and the NS mass observed is taken into account. In our calculations, the 
threshold densities for the appearance of A and U~ are almost the same {pt{A) ~ 
pt{S^) — 4po)j in contrast to the relation pt{U^) < Pt (vl) currently thought to hold. 
This comes from the fact that the repulsive nature of U^n interaction suggested by 
hypernuclear data is duely taken into account (i.e., by using the NHC-D potential). 
Otherwise, we have pt{I^~) < pt{A); the weakening of the U~n repulsion as Vx;-n — *■ 
yj;-„/3, for example, leads to pt{S~) ~ Spo and pt{A) ~ 4.1po- In Table I, we 
list the parameter profiling the y-mixed NSs that correspond to the maximum mass 
-^max and illustrate in Fig. 2 the M-pc relationships indicating how the y-mixed core 
region depends on M and the EOS. We list in Table II the numerical values of the 
composition and EOS for Y-mixed NSs. 

§3. Energy gap equation for hyperon pairing 

As shown in the preceding section, hyperons participate at high densities (p > pt 
(y) ~ 4po), but the fractional density py(= yyp) is relatively low due to small 
contamination, e.g., yy ^ 0.1 for the density region p ~ {pt{Y) ~ 6po) of interest. 
Here, we limit our study to the treatment of hyperon paring in a non-relativistic 
framework. We do this both because the fractional density of hyperons is not high 
and also because, in its present form, the treatment^^^ in the relativistic framework 
involves ambiguities coming from the introduction of a phenomenological cutoff mass 
in the pairing interaction. Corresponding to the small py, the scattering energy EyY 
in the laboratory frame for a {q, -q)-Coopei pair near the Fermi surface {q ~ qpy = 
(Svr^py )^/^ with qpy being the Fermi momentum of Y) is at most 110 MeV, estimated 
by Ep^ = AEpy = Ah'^qj^y/'^My with Epy and My b eing the Fermi kinetic energy 
and the mass of hyperons, respectively. This means that the pairing interaction Vyy 
responsible for y-superfluidity should be that in the ^5*0 pair state {Vyy(^So)) which 
is most attractive at low scattering energies. Thus the gap equation to be treated 
here is of the well-known ^ So-type: 

Ey{q') 

POD 

{q' I Vyyi^So) \ q) = r^drj^{q'r)Vyy{r-^ So)joiqr), (3-4) 
Jo 

including the case of finite temperature with kb being the Boltzman constant. In the 
above expressions, Ay{q) is the energy gap function, €y{q) the single-particle energy 
and My the effective mass of Y. For simplicity, the effective mass approximation 
for ey{q) is adopted in Ea. H3-3|l . We use a bare Vyy{^So) for the pairing interaction 
and do not introduce the G-matrix based effective interaction VVy(^5'o) in place 



^- If ,'%VI ^V-H-5,)l.)f|lt.n.(gM). ,3.1, 
s + AUq'), (3.2) 

= ey{q') - eyiqpy) ^ TTTT^il'^ " Ipy), (3-3) 
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of VyyC^Sq) as done by Balberg and Barnea,'^'^^ because the gap equation plays the 
dual roles of taking account of the pairing correlation and the short-range correlation 
(s.r.c.) simultaneously and hence use of VVy(^S'o) leads to incorrect results with 
larger energy gaps due to the double counting of s.r.c. 

The gap equation fEa. (|3-l() ) is solved numerically with qfy, My and VyyCSq) 
given. We solve it exactly without any approximation, using an iterative technique. 
The QFY (equivalently yy) and My (equivalently the effective mass parameter rriy = 
My /My) are derived in §2. In the following, we make some comments regarding the 
important ingredients in Ea. ()3-1() . my and Vyy{^So), controlling the resulting energy 
gap Ay = AYigpy)- 

i) Effective mass parameter niy 

Generally, the energy gap is very sensitive to my and a larger value is obtained 
for larger my. This is due to the fact that the energy cost for the excitation of 
a yy-pair (ey(g') in Ea. (|3-3|) ) is smaller for larger My(= myMy). In this sense, 
F-superfluidity is more likely to occur than superfluidity of nucleons, since niy ~ 
(0.8 ^ 1.2) is significantly larger than ~ (0.5 ~ 0.7) for p > Apo as shown 
in Fig.3. Moreover, the fact that My ~ (1116 MeV or 1192 MeV)> Mn ^ 940 
MeV strengthens this tendency. Since the single-particle energy has the same effect 
when My = M]^, the effect of niy = 0.8 in the YY pairing corresponds to that of 

= (My /M]\f)my ~ 1.0 in the A^A^ pairing, for the Y = case. A comparison 
of the A and U~ cases reveals that Z'~-superfluidity is more favourable than A- 
superfluidity since 'm*^_ > and in addition M^- > Mj{. That is, the resulting 
energy gaps for baryon superfluids satisfy the relation A^- > Aa > An, as far as 
the role of effective mass is concerned. 

ii) Pairing potential 

The attractive effect of the YY pairing potential Vyy(^S'o) is a main agency for 
the occurrence of y-super fluidity. Our present knowledge of Vyy(^5'o), however, is 
very limited. As in previous works, we adopt three potentials, the ND-Soft^^'^^^'^^) 
**, Ehime^^-* and FG-A,^'^^ with due regard to the present uncertainties of the YY 
potentials. They are commonly based on the one-boson-exchange (OBE) hypothesis 
with SU{3) symmetry in the framework of nonet mesons and octet baryons (B). The 
main differences among them regard the treatments of the short-range interactions 
and the kinds of mesons introduced. The ND-Soft potential is a soft-core version of 
the original Nijmegen hard-core model D (NHC-D)^^) obtained as a superposition of 
three-range Gaussian functions and is constructed so as to fit the t-matrix from NHC- 
D with the hard-core radii being taken as 0.5 fm. The soft-core version is used for 
the sake of convenience in the treatment of the gap equation. The Ehime potential 
is characterized by an application of the OBE scheme throughout all interaction 
ranges and by the addition of a phenomenological neutral scalar meson to take 
the 27r-correlation effects of this sort into account. It consists of a superposition 
of Yukawa-type functions in r-space regularized by the form factors in momentum 
space and has velocity-dependent terms. The FG-A potential is based on the OBE- 

** In the Table I of Ref. 2), V-"^- with i = 2 should be corrected as -2.9232 -> -29.232. 
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Fig. 3. Effective mass parameters m* for baryon components (i = n, p, A, S~) as functions of p 
in the three cases TNI2u, TNI6u and TNI3u for F-mixed NS matter with different stiffnesses. 
The crosses denote the threshold points for the F-mixing. 

scheme, in which the cr-meson is treated in the nonet scheme by taking account of 
the broad width of a- and p-mesons, and it has a Gausian soft-core repulsion with 
strengths constrained by the SU{3) representation of two baryons. This potential 
includes the velocity-dependence and the retardation effects. 

Here we add a comment on the treatment of the channel coupling effect (CCE) of 
AA-EE-SN type in these potentials. The ND-Soft is constructed so as to simulate 
the ^-matrices calculated with the NHC-D in the AA and SN channels. For sim- 
plicity, here, we use only the AA-AA diagonal part of ND-Soft, because the AA-SN 
coupling effect was found to be small for the coupling constant of the model D.*** 
In the Ehime case, CCE is irrelevant, because it is constructed in a single-channel 
approximation. In the FG-A case, however, CCE is significant, and it is taken into 
account by adding an extra term AVsimir) to the direct-channel part V£^(^So), so 

that vjyf\'^So) = Vf^C^So) + AVsimir) simulates the AA phase shifts including 
CCE. In practice, we use this v!{f\^So) as Vaa{^So) for the FG-A potential by 
setting AVsimir) = 93e-(^/'^i)' - lOOOe-^^/'^^)^ ^eV with n = 1.0 fm and r2 = 0.6 
fm. 

*** The statement "this channel coupling effect, although small, is included similarly through the 
original t-matrix to be fitted" appears in Ref. 22) (also similar statements in Refs. 2) and 3)) is 
inaccurate and should be revised as above. 
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Fig. 4. (a) AA scattering phase shifts in degrees as functions of the laboratory frame energy -B^fJi 
in MeV for three AA ^So potentials (ND-Soft, Ehime and FG-A) explained in the text. NN 
scattering phase shifts from the OPEG-^E potential are also shown for comparison. The AA- 
EE- EN channel coupling (CC) effects are indicated by the arrow for the FG-A case (The long 
dash-dotted Uine is for no CC.). (b) Same as (a), but for S~ scattering. 

In Fig.4, the phase shifts from the three potentials are displayed as functions 
of Eyy- For AA, they give similar results at low £^yy, while the discrepancies 
among them increase as EyY increases, which reflects the short-range behavior in 
the YY potentials. For the YY potentials, experimental information exists only for 
the AA case, i.e., double A hypernuclei. It has been confirmed that ND-soft and 
Ehime reproduce the bond energy ABaa{— (4 — 5) MeV) for ^^B.^^)~^^) 
This reproduction is well expected for FG-A because the phase shifts from FG-A 
arc similar to those from ND-Soft and Ehime at low scattering energies. Recently, 
however, there arises the problem that /iByiyi ~ (4 — 5) MeV extracted from ^"^Be 
and ^^B (old data) might be too large and instead ABaa ~ 1 MeV from ^yjHe 
(the "NAGARA event", new data) seems to be correct. This suggests that the 
AA interaction is less attractive than previously thought. If this is true, it has 
a significant implication for the realization of ^-superfluidity. Discussion of this 
problem is given in §5. 

§4. Numerical results and discussion 

4.1. Realization of hyperon superfluidity 

We discuss how the energy gap for hyperon pairing can be realized by referring to 
the well-known case of nucleon pairing. We focus on the effects of pairing interactions 
and effective-mass parameters. For this purpose, the energy gaps Ai (= Ai^qpi); 
i = A and p) calculated from Eqs.(3.1)-(3.4) at T = 0, as functions of the fractional 
density Pi{= Vip), are compared in Fig.5, where the pairing interaction Vii{^So) 
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Fig. 5. Energy gaps Ai for the proton {i = p) 

and Lambda {i = A) as functions of the 
fractional density pi {— yip). This shows 
the effects of three factors, i.e., the pair- 
ing interaction, effective mass M* and frac- 
tional density pi (= yip) controlling the so- 
lutions of the energy gap ctjuation. This 
plot provides an understanding of the real- 
ization of yl-superfluidity. The pairing po- 
tential used for AA (pp) pairing is ND-Soft 
(OPEG ^E-1). 
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Fig. 6. Hyperon energy gap Ai [i = A, 
E~) in y- mixed NS matter, expressed in 
terms of the critical temperature T*^ ~ 
0.57 Ai/ Kb in Eq.(4.1), as functions of the 
total baryon density p. The solid, dashed 
and dotted lines correspond to the ND- 
Soft, Ehime and FG-A potentials, respec- 
tively. The arrow indicates the reduc- 
tion of T^s- when A A interaction is made 
less attractive, as suggested by the "NA- 
GARA event" ^yiHe, by appropriately al- 
tering the baryon-baryon interaction model 
with OBE and SU(3) symmetry frame- 
works. The reduction of Aa is significant 
and the resulting Aa is too small to be de- 
termined, i.e., T*A < 10*^ K. 



is taken from the ND-Soft potential fox i = A and the OPEG ^E-1 potential^''') 
for i = p. We observe the following three points. First, the relation that Ap is 
larger for M* = 0.8Ma than for M* = 0.8Mp means that p-superfluidity is more 
likely if the bare mass of p is as large as that of A. Conversely, yl-superfluidity 
is realized more easily than p-superfluidity due to AIy\ > Mp for the same m*, yi 
and Vii(^So), as mentioned in §3. Second, comparison of Ap with Aa for the same 
effective-mass {M* = M'^ = O.SMa) and the same yi shows that the influence of the 
pairing interaction is to make p-superfluidity more likely than yl-superfiuidity, which 
is reflected by a stronger attraction of Vi\jn{^ Sq), as shown in Fig. 4(a). Third, 
for the same yi, the value of Ap with M* = 0.5Mp is smaller than that of A^ with 
M\ = O.SMyi, in spite of the stronger attraction of V^p(^S'o). This suggests that if the 
actual m*{p) in medium (Fig.3) is taken account (i.e., the case m\ > 0.8 is compared 
to the case m* ~ 0.5 for p > 4po), Ay{ would be much larger than Ap. In fact, as 
shown just below, we have yl-superfluidity and also Z"" -superfluidity in NS cores, 
but we have no ^^-superfluidity due to a significant effect coming from the very small 
m*{p) and relatively large pp (namely yp). Thus, we can say that y-superfiuidity in 
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dense NS cores occurs as a result of the fact that the attractive effect of Vyy{^So) 
is not so different from that of VNNi^So) and is mainly due to the strong influence 
of my(p) which is much larger than mj^{p). 

Results for a realistic case with p-dependent my and yy are displayed in Fig. 6 
for the TNI6u NS model as an example, in terms of the critical temperature T*y 
for y-superfiuidity defined by the well-known formula from the near-Fermi-surface 
approximation, 

T*Y = 0.57 Ay/ KB ^ OMAy x W^°K, (4-1) 

with Ay in MeV calculated at T = {Ay = Ay{qF, T = 0)). The solid, dashed and 
dotted lines correspond to IVy(^S'o) from the ND-Soft, the Ehimc and the FG-A 
potentials, respectively. The y-superfluidity occurs when Tcy exceeds Tj„ ~ 10^ K, 
the internal temperature of usual middle-aged NSs. The following points are noted: 



Table III. Energy gaps Ay of A and in TNI6u y-mixed NS matter at zero temperature {E = 0) 
for several densities p and for three YY pairing interactions; ND-Soft, Ehime and FG-A. The 
mixing ratios yy and the effective mass parameters my(= My /My) of hyperons in medium are 
also given. Ay at finite temperature (T > 0) is obtained from the profile function P(r) given 
in S4.2. 



Y 


p/po 


VY 


my 


A 


y in MeV 










ND-Soft 


Ehimo 


FG-A 




4.05 


0.00057 


0.818 


0.006 


0.019 


0.006 




4.10 


0.00225 


0.820 


0.061 


0.133 


0.067 




4.15 


0.0046 


0.821 


0.140 


0.274 


0.168 




4.25 


0.0102 


0.825 


0.274 


0.496 


0.370 




4.50 


0.0261 


0.832 


0.339 


0.638 


0.635 


A 


4.75 


0.0422 


0.839 


0.219 


0.486 


0.641 




5.00 


0.0578 


0.845 


0.089 


0.271 


0.516 




5.10 


0.0638 


0.847 


0.053 


0.196 


0.450 




5.25 


0.0726 


0.850 


0.019 


0.106 


0.348 




5.35 


0.0782 


0.852 


0.008 


0.065 


0.285 




5.50 


0.0865 


0.854 


0.001 


0.025 


0.195 




6.00 


0.1116 


0.860 


0. 


0. 


0.027 




4.10 


0.00156 


1.165 


1.656 


2.295 


15.866 




4.15 


0.0047 


1.146 


2.428 


3.166 


16.502 




4.25 


0.0123 


1.121 


3.238 


4.125 


17.442 




4.50 


0.0321 


1.087 


3.619 


4.778 


18.440 




4.75 


0.0506 


1.062 


3.127 


4.514 


18.423 




5.00 


0.0677 


1.040 


2.300 


3.867 


17.901 




5.25 


0.0833 


1.019 


1.385 


3.037 


17.062 




5.50 


0.0975 


1.000 


0.612 


2.165 


16.036 




5.75 


0.1104 


0.983 


0.154 


1.362 


14.895 




6.00 


0.1222 


0.967 


0.011 


0.714 


13.649 




6.50 


0.1425 


0.938 


0. 


0.077 


11.010 




7.00 


0.1591 


0.914 


0. 


0.000 


8.326 
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Fig. 7. Hyperon energy gap Ay for the ND-Soft pairing potential as a function of p for TNI2u, 
TNI6u and TNI3u EOSs. This reveals the strength of the EOS dependence of Ay. The left-hand 
(right-hand) side depicts the A{E~) case. 



(i) Both A- and Z'^-superfluids are likely to exist in NS cores, since TcA and T^x:- 
are well above Tj„. A and U~ become superfluids as soon as they appear in NS 
cores, and yl-superfluidity is realized in a limited density region (/? ~ (4 — Q)po), 
while i7~-superfluid extends to much higher densities. 

(ii) The behavior of the Tcy versus p relationships for three different potentials are 
more closely bunched in the A case than in the U~ case, which comes from the fact 
that there is less uncertainty in the pairing interaction in the A case, due to the 
experimental information from double A hypernuclei mentioned in §3. 

(iii) T^x- (~ 10i°-ii K) for i;~-superfluid is larger than T^a (~ W^-^^ K) for 
yl-superfluid by more than one order of magnitude, which comes mainly from the 
very large (~ 1) compared to m\ (~ 0.8). The large T^^- implies that U~- 
superfluidity can be realized in NSs even at an early stage of the thermal evolution, 
when Tin is as high as 10^° K. 

The full results for the hyperon energy gaps for TNI6u NS models and for 
three VyyCSq) are presented in Table III at several densities, together with the 
p-dependent input parameters, yy and my, controlling the solution of the energy 
gap equations H3-l|) - H3-4|) . In the cases of TNI3u and TNI2u, the points (i) - (iii) 
noted above are very similar to those in the TNI6u case. Thus, we can conclude 
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that A and admixed in NS cores are certainly in the superfluid state, although 
this depends somewhat on the NS model and considerably on the pairing interaction 
VyyCSq). The EOS-dependence of hyperon gap is shown in Fig.7. 

4.2. Effects of finite temperature 

So far we have investigated the energy gaps at zero temperature. In cooling 
calculations of NSs, however, we need to take account of the temperture effects 
on the energy gaps, as the temperature T of Y-mixed NS matter varies during 
the thermal evolution and in general the resulting energy gap Ay decreases with 
increasing T.^^^ The T-dependent energy gaps Z\y(T), the solution of the finite- 
temperature gap equations p-ip - H3-4() . are plotted in Fig.8 for a typical case with 
ND-Soft and TNI6u. We observe that Ay depends remarkably on T, especially for 
small Ay. In the cooling calculations, it is of practical use to have a profile function 
P{t) approximating the T-dependence of Ay. Following our previous work,^^) we 
introduce P{t) as a function of the argument r defined by r = T/Tc'. 



Py{T) = Ay{qFY,T)/Ay{qFY,T = 0), (4-2) 

= 1 - ^/2^i^exp[-l/(aor)](l + air + a2r2), {0 < t < n) (4-3) 

= aoCi^/r^[l + /3i(l-r)+/32(l-r)2], {n < t < 1) (4-4) 

ao = a*{TjT:), (4-5) 

where oq, Ci, ai, 02, Pi and P2 are constant parameters and Oq = 0.57 is given by 
the coefficient appearing in the definition of T* in Eq. (|4Tj) . Here Tc (T* ) denotes the 
actual (approximate) critical temperature. Equation H4-3() ensures that the energy 
gap ratio (Ea. H4-2() ') decreases from P{t) = 1 as T starts from T = 0, and Ea. H4-4() 
ensures that P{t) = at T = Tc. 

From numerical calculations used to search for Tc, it is found that Tc is well 
substituted by T*, except in the case of S~ with the FG-A potential, where, due 
to the particularly large Ajj-, we have Tc ~ (0.278 + 0.130{p/ po)) x T* at densities 
(4 — 6)po- By selecting the boundary value of the two regions of as Tb = 0.7 
and fitting A{qpy, T = O)P(r) to A{qFy, T) actually calculated, we obtain the 
parameters of P{t) listed in Table IV (details of the procedure are given in Ref. 49)). 
The reproduction of A{qFY-, T)/ A{qFy, T = 0) by P{t) is illustrated in Fig. 9 for 
A as an example. Once P{t) is obtained, we can use an analytic form of the T- 
dependent energy gap in NS cooling calculations by utilizing the zero temperature 
gap Ay{qFY, T = 0). 



Table IV. Parameters for the profile function P(r) (= Ay{T)/Ay{0)) analytically representing the 
temperature dependence of the hyperon(Y) energy gap. 



Y 


Ci 


ai 


a2 


Pi 


132 


A 


2.8963 


-0.45902 


1.41171 


-0.13926 


-0.53396 


E- 


2.8504 


-0.39488 


1.38328 


-0.09447 


-0.57033 
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Fig. 8. Temperature (T) dependence of the A 
energy gap (Aa) for the case of the TNI6u 
EOS and the ND-Soft potential, as an ex- 
ample. 




Fig. 9. Profile fimction Pa (t) (= Aa{T)/Aa{0)) 
in the form of Eas. l|4-8ll and II4-4II as a func- 
tion of T (= T/Tca), averaged over the pair- 
ing potentials employed and the densities 
under consideration (p ~ (4 — 5)po)- The 
crosses and the error bars denote the val- 
ues obtained from numerical calculations of 
the finite-temperature gap equation. Here, 
Tt — 0.7 denotes the boundary between the 
region of Eq. (g^ and that of Eq. . The 
parameters for Py(t) {Y = A and E~) are 
listed in Table IV. 



4.3. Momentum triangle condition for Y-DUrca cooling 

The z^-emission process responsible for NS cooling acts only when a certain con- 
dition, namely, the so-called "momentum triangle condition" relevant to momentum 
conservation in the corresponding reaction, is satisfied. For example, a familiar (3- 
decay process associated with nucleons (A^) and its inverse process (called A^-DUrca), 
i.e., p^n + l + Vi^n + l^p + vi (abbreviated as p <-> n/ hereafter), is allowed 
when the Fermi momenta ki = qFi{= (Svr^yj/?)^/^) for particles participating in the 
reaction satisfy the triangle condition 

\ kp - ki \< kn < kp + ki. (4-6) 

Here the Fermi momentum ki, (~ K-QT/hc) for neutrinos is omitted, because it is 
much smaller than those of N and I. As is well known, A^-DUrca is forbidden in 
usual neutron star matter including n, p, e~ and /i^ components, because the sec- 
ond inequality in Ea. ()4-6() is hard to hold among n (large component, y„ > 0.9), 
p (small component, yp < 0.1) and / (small component, yi < 0.1). Then, to sat- 
isfy the condition, it is necessary for another particle to participate as a bystander 
mediating momentum conservation. Indeed, z^-emission due to the modified Urea 
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process {pN <-> nNl) including a bystander nucleon becomes a main agency for the 
cooling of usual NSs, although the z/-emissivity is much smaller than that of the N- 
DUrca. In the case of normal NS matter {npe~ fj,~ matter), the fast cooling due to 
A^-DUrca is made possible only for NS models including high proton fraction phase 
with yp > 0.148.5°) 




Fig. 10. Momentum triangle condition to check the functioning of various DUrca i/-emission pro- 
cesses discussed in the text, where ki — qpi (= (37r^t/ip)^^^) is the Fermi momentum of the 
respective constituents {i = n, p, A, e~ , n~) in F-mixed NS matter with TNI6u EOS; the upper 
panel depicts the n *-* p + 1 process, the middle panel the S~ ^ n + I process, and the lower 
panel the S~ yl + / process, depending on p. The solid (dotted) lines corresponds to the case 
with / = e~ {n~ ) . 

How about the situation in the y-mixed NS matter {npAS~e''ii~ matter) under 
consideration? In Fig. 10, the triangle condition is checked for DUrca processes, i.e., 
A^-DUrca and Y-DUrca, for the K-mixed NS matter of TNI6u EOS by using the 
fractions yi listed in Table 11(b). The behavior for other EOS, the TNI2u EOS 
(Table 11(a)) and TNI3u EOS (Table 11(c)), are very similar to those in Fig.lO. We 
note the following points: 

(i) In the F-mixed phase, the proton fraction yp becomes remarkablely larger than 
that of usual NS matter {npe~ ji~ matter), increasing with p and reaching 20% at 
p ~ 8/9o. However, Fig.lO shows that A^-DUrca does not function up to the high 
density region, i.e., p < 6.5po for n <-> pe~ and p < 8po for n ^ PfJ'~- Combining 
this with the M-pc relationship in Fig.2, this means that AT-DUrca cooling begins 
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to function only for rather massive NSs with M > 1.6Mq for n ^ pe DUrca and 
M > 1.7 Mq for n ^ pii~ DUrca. 

(ii) In contrast with the A'^-DUrca mentioned above, the y-DUrca of yl pZ type 
can operate as soon as A begins to appear, i.e., for p > 4po both for / = e~ and 
I = li~ . This is because the triangle condition can be satisfied readily due to the 
fact that all the particles participating in the reactions belong to small components, 
not including the large n-component. From Fig. 2, we see that A ^ pi y-DUrca can 
work in all the F-mixed NSs with M > 1.35Mq, and it provides a very fast cooling 
of NSs. 

(iii) It is difficult to realize conditions under which the y-DUrca of E~ nl type 
including the large n-component functions. In fact, this is impossible for NSs with 
central density pc < 7po for E~ ^ ne~ and p^ < S.Spo for ^ np~. This implies 
that IJ~ •i-^ nl y-DUrca occurs only for NS with M > I.GMq in the / = case and 
M > 1.7Mq in the I = fi~ . On the other hand, like the A ^ pi type, the y-DUrca of 
the U~ ^ Al type occurs as soon as E~ appears, because all the associated particles 
are small components. It should be noted, however, that the y-DUrca including E~ 
is made possible only for the density region where A coexists with E~ . That is, the 
participation of A in NS constituents is essential for the y-DUrca to function. 

Based on (i) - (iii), it is remarked that the momentum triangle condition inter- 
estingly limits the density region (in other words, NS mass) and the type of DUrca 
processes. Except in the case of massive NS stars, A mixing is a main mechanism 
for y-DUrca cooling. 

4.4. Superfluid suppression effects on neutrino-emissivities 

As mentioned in §1, the fast cooling due to y-DUrca alone is not sufficient 

to explain colder class NSs, because its extremely efficient i^-cmission causes the 
serious problem of too rapid cooling. Here we discuss the suppression effects on v- 
emissivities due to nucleon and hyperon superfluidities, using the energy gaps in 
Table II and their T-dependcnce given by the profile function P{t) in §4.2. The 
y-DUrca i/-cmissivities, eD\j{A ^ p) = eDu(^ ^ P^ ) + (A ^ PIJ'~) and eDu(^~ 
A) = eDu(^~ ^ Ae~) + eDu(^~ ^ ^A*~)) given by^^^ 

€Du(^ ^p) = 0.158 X lO'^TKyeP/po)'/^ 

xm\m;{MAMp/M^)RAA{A,p)i9e + 9^), (4-7) 

eDu(^" ^ .1) = 0.822 X lO^irKi/ep/po)'/' 

xm*^^m\{Ms-MA/Ml)RAA{S-,A){ee + 6^), (4-8) 

where Tg = T/W^ K, 9i = 1 (0) for the process allowed (not allowed) and Raa{Bi, 
B2) is the superfluid suppression factor relevant to the participating baryons Bx and 
B2. The appearance of y^, comes from the chemical equilibrium; /Xe = /"^ ~ hqp(.c. 
The quantity Raa{Bi, B2) is defined by the ratio of the integral of the statistical 
factors in the ^Sq superfluid phase (Ig) to that in the normal phase (/q), i.e., Raa{Bi, 
B2) = Is/Iq, which represents the decrease of the number density around the Fermi 
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surface due to the energy gap (further details including calculational methods are 
given in Ref. 49)). 



Log 8 (erg/cm^s) 

LoqRaa 26 1 ' ■ ■ , ^ 




5 p/po 6 



Fig. 11. (a) Superfluid suppression factor Raa{Bi, B2) for !/-emissivity relevant to the baryons (Si, 
B2) participating in the F-DUrca processes, as functions of p and for two typical temperatures 
of the NS interior, T = 1 x 10* K and 5 x 10® K. The sohd lines (dashed lines) is for A <-> p. 
{E~ ^ A) process summed over A ^ p + er and A ^ p + l-i^ {S^ ^ A + e" and ^ A + ij,~). 
(b) y-DUrca i^-emissibity including the superfluid suppression effect corresponding to Raa{Bi, 
B2) in (a), as functions of p and two typical temperatures. The solid lines (dashed lines) is for 
y-DUrca of A *~* p {S~ ^ A) type and the dotted lines are for the modified Urea processes 
responsible for the standard cooling of NSs. The emissivitics without the superfluid suppression 
are also included as a reference and are plotted by the dash-dotted lines in the case of the A p 
process. Calculations in (a) and (b) are performed for TNI6u EOS and energy gaps from the 
ND-Soft pairing potential. 

Calculated results are shown in Fig. 11 for TNI6u EOS and energy gaps from 
the ND-Soft pairing potential. Figure 11(a) demonstrates the suppression effects 
Raa{Bi, B2) as a function of p by considering two typical NS temperatures, T = 
1 X 10^ K and 5 x 10^ K. Raa{.Bi, B2) depends strongly on both p and T, as reflected 
by the p- and T-dependent energy gaps of Bi and B2- For example, Raa{-^i p) for 
A p process (solid lines) depends on p as Raa ~ 1 ^ 10^^ 1 for p/po ~ 
4 — > 4.5 — 5.0 at r = 5 X 10^ K. When the temperature decreases to T = 1 x 10^ 
K, the suppression is much stronger; Raa — 1 ^ 10^^^'^ 1. Apart from the 
details, this characteristic feature can be qualitatively understood as resulting from 
an exponential damping factor, exp[— (Z\bj + AB2) / hbT] involved in Raa{Bi, B2), 
together with the p- and T- dependences of Ab^ and in Fig. 6 {in A p process, 
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Ap = as mentioned in §4.1.). For the S <-> yl process, the p- and T-dependent 
behavior of Raa{^~ , ^) is extremely enlarged because of the large energy gaps of 

Figure 11(b) illustrates how the emissivities eDu(^i ^ B2) for y-DUrca are 
moderated by Raa{Bi, B2), where the emissivities eMD for the modified Urea in- 
cluding superfluid suppression (see Ref. 49) for details) are also shown for compari- 
son. First, wc consider e{A p) (solid lines) in the density range p ~ (4 — 5)po- At 
T = 5 X 10^ K, the effect of Raa{j^, p) is not significant and eDu(^ ^ p) ~ 10^3-24 
erg/cm^s is much larger than eMU ~ loi^"^^ erg/cm^s. However, at T = 1 x 10^ 
K, the situation is very different; eDu(^ p) ~ 10^° erg/cm^s for no suppression 
(dash-dotted lines) is maximumly suppressed to eDu(^ ^ p) ^ 10^'^ erg/cm^s by 
Raa{-^-, p) and is even smaller than cmu ~ 10^^^^"^ crg/cm^s for p ~ (4.2 — 4.7)/3o, 
although €du(^ ^ p) ^ ^MU continues to hold for other densities. The behavior of 
eDu(^ p) for T = 5 X 10^ K — 1 X 10^ K reveals that at higher T in relatively 
earlier stage of the thermal evolution, the A ^ p process dramatically accelerates 
the NS cooling as compared to the standard modified Urea process, but this accel- 
eration experiences a braking effect due to superfluid suppression as T decreases in 
later stage, and in this way, too rapid cooling does not occur. In the case of the 
E~ ^ A process, eDu('^~ ^ ^) is much smaller than emU) as well as eDu(^ ^ p)i 
for p ~ (4 — 5)p0) due to the extremely strong suppression caused by Raa{^~ , A). In 
the total z^-emissivity, the most efficient i/-emission process is of prime importance. 
So the ^ A process becomes appreciable only at high densities (p > 5.5po) where 
eDu(-^^ ^ A) becomes comparable to or larger than eDu(^ ^ p)i depending on T. 
At high densities (p > 5po)) eDu(^ ^ p) and/or eDu(-^~ ^ ^) suffer no suppres- 
sion, because of the disappearance of the superfluidities, and hence massive NSs 
with pc > 5po cause a problem of too rapid cooling. This undesirable situation is 
the same as that for the A^-DUrca cooling mentioned in §4.1 which becomes possible 
for p > 6.5po, because the n and p superfluidities completely disappear for such high 
densities. 

Taking notice of the two conditions, i.e., the triangle condition to allow fast 
DUrca cooling and the superfluid suppression condition to moderate the cooling 
rate (avoiding too rapid cooling), we can say that the y-DUrca of yl <-> p type 
plays a decisive role in bringing about the cooling scenario for colder class NSs.^^) 
This is because firstly the ^ A process can function only with the coexistence 
of A and secondry the large gap suppresses e-£)\j{E~ <-> A) too strongly . To 
obtain consistency with surface temperature observations of NSs, we can constrain 
the mass of the colder class NSs as follows. For the TNI6u NS models, it must be 
larger than about l.SSM© in order to realize pc > Pt(X) and also be smaller than 
about I.55M0 in order to avoid too rapid cooling. Of course this mass restriction 
depends on the NS model (Fig.2) and the superfluid gap A (Fig.6), but in any case 
the functioning of the y-DUrca moderated by the baryon superfluidity provides us 
with a new idea regarding NS cooling. That is, less massive NSs (hence without 
y-mixed core) cooled slowly by the standard process (modified Urea) correspond 
to hotter class NSs, and massive NSs (with y-mixed core and superfluids) cooled 
rapidly by the nonstandard process (y-DUrca) correspond to colder class NSs. In 
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fact, cooling calculations (although of a preliminary nature), nicely account for the 
surface temperature observations employing this idea.^^^ 

§5. Less attractive AA interaction 

We have discussed that the hyperon cooling scenario for NSs is consistent with 
observations and the occurrence of Tl-superfluidity is of particular importance. There 
exists, however, an experimental result that in incompatible with this occurrence. 
Recently the KEK-E373 group observed one event of the A hypernucleus ^^jHe in an 
emulsion experiment, called the "NAGARA event" and extracted new information 
of ABaa — 1 MeV for the AA bond energy from the difference between the binding 
energies B of double and single A hypernuclei:^^^ 

ABaa = 5(%He) - 2 x B{\Re). (5-1) 

The new result ABaa — 1 MeV is much smaller than the old one ABaa — (4 — 5) 
MeV extracted from ^*^Be and ]\a^- This suggests that the AA attraction responsible 
for yl-superfluidity would be much weaker than previously thought. Since all the 
pairing inter acations, ND-Soft, Ehime and FG-A, used in the present calculations of 
the ^-gap (§3) are taken to be consistent with the old ABaa — (4 — 5) MeV, the 
new finding of ABaa — 1 MeV must be considered seriously and it necessitates a 
reexamination of yl-gaps for the reduced AA attraction. Hiyama et al.,^^-* performing 
energy calculations for the a+yl+yl system, found that ABaa — 1 MeV is reproduced 
by taking V^^{^So) V^^{^So) ~ 0.5V^J^{^So). More precisely, the potential 
strength in the modified ND-Soft y^(^S'o) is multiplied by a factor of 0.45 for 
the shortest-range term and a factor of 0.5 for the intermediate- and long-range 
terms. By using this V^(^S'o), we have tried to calculate A a and found that yl- 
superfluidity does not occur, i.e., Tc(A) << Tj„ ~ 10^ K. This leads us to inquire 
about the situation for Ajj- since YY interactions are internally related according 
to the OBE plus SU{3) symmetry hypothesis. To see this, we reconstruct the BB 
interaction model of FG-A so as to reproduce the AA phase shifts from V'^(^'S'o) 
by weakening the attractive contribution form the cr-meson exchange. Through new 
parameters adjusted, new U~ U~ interaction V^-j^- is to be obtained. This can be 
done as follows. 

In the framework of octet baryons plus nonet mesons with SU (3) symmetry, the 
coupling constants (jBBm = dBBm/V^ for the BB interaction by single meson (m) 
exchange are given in terms of four parameters, i.e., the singlet coupling g^^\ the 
octet coupling g'^^^ , a parameter a relevant to the so-called F-D ratio, and the mixing 
angle 6. Therefore, the set of these four parameters is uniquely determined when 
four coupling constants are given. Here we stipulate that the NN interaction sector 
be unaltered as it is firmly established. There we use, as four coupling constants, 
dNNcr, dNNao 9NNfo which are unchanged from the original FG-A model, and 
9AAa which is adjusted so as to reproduce the phase shifts from V^(-^S'o) compatible 
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Fig. 12. A demonstration of to show how the attractive features of ^So YY interactions (Fig.4) are 

weakened by taking account of the less attractive AA interaction suggested by the "NAGARA 
event" ^yiHe in the framework of OBE plus SU{3) symmtry. The solid lines (dashed lines) are 
the phase shifts for the original FG-A potential Vyy (the modified FG-A potential Vyy)- 
The modification {Vyy ^ Vyy) has been done through the fitting of Vj^ to VaT obtained by 
seeking reproduction of the AA bond energy of aaH-^- The phase shifts from Vaa sue shown by 
the dotted line for reference. 

with the experimental results for ^^He: 

^jviva = g^^^ COS 6 + {4a- 1) /Vs ■ g^^^ sin 6, (5-2a) 

gNNao = g^^^ (5-2b) 

gNNfo = (4a- l)/\/3-5(^)cos6>-5(^)sin6', (5-2c) 
gAAa = g^^^ COS 61 - 2(1 - a) /V3 • 5^^^ sin 6. (5-2d) 

Actually, about 10% reduction of cjAAo- is satisfactory for the reproduction of 
the phase shifts at low scattering energies ( < 50 MeV) corresponding to the pairing 
problem of A. Prom the new set {g^^\ g^^\ a, 9}, new gAAm and (jE-s-m responsible 
for the modified VJ^Q-Sq) and y|'^^_(^S'o) are obtained as in Table V. Figure 12 
compares the attractive effects of V^yy (^5*0) (solid lines) and Vyy(^So) (dashed lines) 
in terms of their phase shifts. We see that the S~ E~ attraction of modified FG-A, 
as in the AA case, is significantly reduced from that of the original FG-A. However, 
V^'!fj,_Q-So) remains to be more attractive than the original V^(^S'o), implying the 
persistence of ^'"-superfluidity. In fact, a reexamination of A^- with V^9^_(^S'o) 
shows that i7~-superfiuidity is realized with Tc{I^~) reduced by about a half, as 
indicated by an arrow in Fig.6. Also note that Fig. 12 shows that the situation for 
the E~E~ case is similar to the E~E~ case, suggesting that E~ when admixed in 
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NS cores could be in a superfluid state. 

To summarize, the less attractive A A interaction suggested by the NAGARA 
event leads to the disappearance of yl-superfluid, although X!~ superfluidity can per- 
sist, and the y-cooling scenario encounters the serious problem of too rapid cooling 
because the superfluid suppression of the i^-emissivity for the A pi process does 
not occur. 

Table V. Coupling constant gvYm for the modified FG-A potential {Vyy ) to take account of "NA- 
GARA event" data. The numbers in the parentheses are for original FG-A potential {Vyy)- 



gYYm\m 


a 


/o 


Oo 


9AAm 


4.656 


2.0849 







(5.15543) 


(2.46243) 


(0) 




4.4265 


-2.5360 


5.5307 




(5.53592) 


(-1.41267) 


(4.89611) 


gS-S-m 


4.3610 


-3.8542 


4.7687 




(5.66490) 


(-2.72596) 


(4.13408) 



§6. Summary and remarks 

We have calculated the energy gaps of A and Z"" admixed in NS cores for val- 
ues of the mixing ratios and effective mass parameters of the F-mixed NS models 
obtained realistically using G-matrix-based effective interaction approach, and by 
adopting several pairing potentials to take account of the uncertainties on the YY 
interactions. It is found that both A- and i7~-superfiuids are realized as soon as they 
begin to appear at around 4y9o, although the critical temperature and maximum den- 
sity for the existence of the gaps depend considerably on the pairing potential. The 
EOSs for y-mixcd NS matter with different stiffnesses, together with the fractions 
of the respective components, and also the A- and Z'^-energy gaps were numerically 
determined and tabulated. The profile function P(r) is presented to give analyti- 
cally the temperature dependence of energy gaps. The presentation of these physical 
quantities serves as physical inputs for the calculations of Y- cooling for NSs. 

Occurrence of Y superfluidity is indispensable to the y-cooling scenario to be 
consistent with observations of colder class NSs, because otherwise enormous v- 
emissivities due to y-DUrca processes cannot be moderated by superfluid suppres- 
sion effects and this would cause a serious problem of too rapid cooling generally 
encountered in DUrca processes. The momentum triangle conditions to discriminate 
the functioning DUrca processes were checked in y-mixed NS matter. By combining 
the effects of superfluid suppression, it was found that in the y-DUrca processes only 
the A ^ pi and E~ •s-^ Al processes can function, and the former plays a decicive role 
in y-cooling, whereas the latter is completely supressed due to the large Z!~ gaps 
except in the particulary dense core of NSs. The AT-DUrca in its currently accepted 
form is excluded for the fast cooling scenario, because it does not take place in core 
densities of medium mass NSs and its functioning in massive NSs encounters the 
problem of too rapid cooling owing to the complete disappearance of nucleon en- 
ergy gaps. On the basis of the investigations presented here, we conclude that there 
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is a strong possibility that the F-coohng scenario combined with y-superfluidity 
is a promising candidate for nonstandard fast coohng of NSs. This scenario can 
explain observations according to which the hotter class NSs have lighter masses 
(e.g., M < 1.35Mq for TNI6u plus ND-Soft and no K-mixcd core) and are cooled 
slowly by the standard modified Urea, while colder class NSs have medium masses 
(M ~ (1.4— 1.5)Mo with F-superfluids) and are cooled rapidly by a moderated non- 
standered cooling of y-DUrca. The massive NSs (M > 1.55Mq) including a dense 
core with the disappearance of y-superfluids causes too rapid cooling and are ex- 
cluded if extremely cooled NSs are not found. It is worth noting that observations 
of the surface temperature interestingly limit the mass of NSs. 

The less attractive AA interaction inferred from the NAGARA event ^^He affects 
so much the y-cooling scenario. It was shown that yl-supcrfluidity disappears when 
this weaker attraction is taken into account although Z'^-supcrfluidity survives. This 
leads to the consideration that the y-cooling scenario, where yl-superfluid plays an 
essential role, breaks down due to too rapid cooling. However, we cannot draw a 
final conclusion in this regard until the less attractive AA interaction not to allow 
yl-superfiuidity is confirmed by more extensive studies on the following issues: the 
validity of extracting the AA bond energy of a ^yiHe system using Eq.(5.1) without 
considering the rearrangement effects, the possibility of a small ABaa as a result 
of the three-body force repulsion, the mass-number dependence of ABaa in double 
A hypernuclei, possible effects of AA-EE-EN channel coupling in the AA pairing, 
yl-superfiuidity due to AE^ pairing not restricted to AA pairing, and of course the 
new finding of double A hypernuclei. For example, it is of importance to note that 
the inclusion of rearrangement effects on the a-particle in the a + A + A system 
suggests a stronger A A attraction than V^(^5o).^^^'^^^ 

Finally, it is worth while mentioning another possibility for the fast cooling sce- 
nario. As is well known, several processes of nonstandard fast cooling of the DUrca 
type have been proposed. These function under the possible presence of nonstandard 
components, such as a pion condensate ( pion cooling) , kaon condensate (kaon cool- 
ing) , high proton- fraction phase ( A^-DUrca cooling) and quark phase (quark cooling) . 
Among these, A-DUrca and kaon coolings are excluded from the candidates for fast 
cooling scenario consistent with observations, because baryon superfluidity is un- 
likely in these exotic phases.^^^ Also quark cooling is excluded because of too large 
energy gaps exceeding several tens of MeV which completely supress the i^-cmission. 
Then, only pion cooling remains as another promising candidate because there is a 
strong possibility of quasibaryon superfiuids due mainly to the large effective masses 
and pairing attraction enhanced by isobar z:\(1232) contamination.^^)' 
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